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In deposition with a poisoning species, we show that the transition to a blocked or pinned phase 
may be viewed as an absorbing transition in the directed percolation (DP) class. We consider a 
ballistic-like deposition model with an active and an inactive species that represents its basic features 
and exhibits a transition from a growing phase to a blocked or pinned phase, with the deposition 
rate as the order parameter. In the growing phase, the interface width shows a crossover from 
the critical W ~ t behavior to Kardar-Parisi-Zhang (KPZ) scaling, which involves DP and KPZ 
exponents in the saturation regime. In the pinned phase, the maximum heights and widths scale as 
H s ~ W a ~ (p — Pc)~ Vli ■ The robustness of the DP class suggests investigations in real systems. 

PACS numbers: 05.50.+q, 64.60.Ht, 68.35.Ct, 68.55.Ln, 81.15.Aa 



During some deposition processes, the presence of dif- 
ferent chemical species improves films' properties but 
may also lead to undesired features, such as the decrease 
of growth rates due to erosion processes or the saturation 
of dangling bonds at the surface. One important exam- 
ple is the deposition of Si films doped with P by CVD 
or MBE in atmospheres with phosphine Jl]-|j|], in which 
it is observed the decrease of growth rates when phos- 
phine flux increases. This feature seems to be related 
to the saturation of dangling bonds at the surface ||. 
Similar poisoning effect appears in diamond CVD in at- 
mospheres with boron and nitrogen Q. High fluxes of 
the poisoning species may cancel out the growth of the 
main species, thus showing a transition from a growing 
phase to a blocked or pinned phase. Here we will ar- 
gue that, in the absence of erosion processes of these two 
species, it may be viewed as a transition to an absorbing 
state in the directed percolation (DP) class |3|-|| , and we 
will present a deposition model that represents the main 
features of this process. 

We will consider a statistical model that represents the 
essential aspects of films' growth and may be used to 
calculate growth rates, analyse surface roughness scaling 
and predict a dynamic transition. It is a ballistic-like de- 
position model with two species, an active one (A) and 
an inactive one (B), with a continuous transition from a 
growth phase to a blocked phase. The mapping of this 
transition onto the DP class shows that the growth ve- 
locity is the order parameter of the problem and that 
the growth phase corresponds to the active phase of DP. 
Thus their physical properties are completely different 
from previous models of surface growth with pinning or 
roughening transitions J9|-[l3| . The observed fall of depo- 
sition rates in the growth regime agrees qualitatively with 
deposition experiments showing poisoning effects. Thus, 
the interpretation of the pinning process as a transition 
to an absorbing state and the robustness of the DP class 
strongly suggest that other transitions to blocked phases 
due to poisoning of films' growth are also in the DP class. 



Furthermore, we will show that the scaling of quantities 
such as growth rates, surface roughness and thicknesses 
of blocked deposits involve the exponents of the Kardar- 
Parisi-Zhang (KPZ) theory [|l4| and DP exponents, and 
may eventually be used to compare our theory with ex- 
perimental data. 

In the following we will describe our model, show the 
results in one-dimensional substrates while discussing the 
relation to DP, show some results in two dimensions and 
present a final discussion. 

In our model, particles A and B are released from ran- 
dom positions above a d-dimensional surface of length L 
with probabilities 1 — p and p, respectively. The incident 
particle follows a straight vertical trajectory towards the 
surface. Aggregation is allowed only if the incident par- 
ticle encounters a particle A at the top of the column of 
incidence or at the top of a higher neighboring column. 
Otherwise the aggregation attempt is rejected. Fig. la 
illustrates the aggregation rules. A column in which ag- 
gregation is possible will be called an active column. The 
deposition time is the number of deposition attempts per 
substrate column, thus the deposition rate (number of de- 
posited particles per unit time) is equal to the fraction 
of active columns. 

It is clear from Fig. la that particles B represent im- 
purities that prevent the growth to occur in their neigh- 
borhoods. This model resembles the AC model proposed 
by other authors [|l5|jl(| , but their results are very differ- 
ent from ours (a morphological transition was suggested 
in d = 2 ]l6| ] , but it was not quantitatively studied) . Our 
findings are also completely different from the two-species 
RSOS model of Ref. Jl?]], although the pure case (p = 0) 
also obeyed KPZ scaling. 

Now we will present results in d = 1. 

For small values of p, the growth process continues 
indefinitely, such as in the pure model (p — 0). However, 
when p increases, the growth rate r decreases due to the 
increase in the density of B at the surface, as shown in 
Fig. 2a. In Fig. 2b we show lnr versus ln(p c — p) for 



1 



p c = 0.20715, which gives the best linear fit of the data 
for 0.19 < p < 0.206. Thus we obtain 



e^e =Pc 



P, 



(1) 



with p c = 0.20715 ± 0.00010 and f3 = 0.282 ± 0.012. 

The instantaneous growth rate decays as the density of 
particles A at the surface. Focusing on the surface config- 
uration, we notice that the growth rules of Fig. la may 
be mapped onto a d-dimensional contact process [|18||7| 
(CP) in which a top A represents a particle and a top 
B represents a hole (or empty site), as shown in Fig. 
lb. When the deposition of a B occurs in a column with 
a top A, it corresponds to the annihilation of a parti- 
cle in the CP. On the other hand, the deposition of an 
A in a column with a top B and a neighboring column 
with a top A corresponds to offspring production in the 
CP. Notice that the stability of the absorbing state is 
represented by process 4 in Fig. lb. The probabilities 
of annihilation and offspring production in the CP are 
not trivially related to p, since they also depend on the 
neighboging heights' distribution. 

The equivalence to a CP indicates that the transi- 
tion is in the DP class, with the density of top A or 
the growth rate r as the order parameter. The above 
value of the exponent (3 and forward results support 
this statement (the best known estimate for DP is (3 = 



0.276486 ± 0.000008 |19|) 



Here it is relevant to recall that all known statisti- 
cal models showing continuous transitions to absorbing 
states, with positive one-component order parameters, 
short-range interactions and no additional symmetries, 
are in the DP class P,p0|. This so called robustness of 
the DP class is the reason for us to expect universal- 
ity in real systems' transitions with the same blocking 
mechanisms of our model. 

At the critical point in d = 1 and L < 8192, we esti- 
mated the deposition rate at very long times, r^L), and 
obtained 



. (£, t = oo ) ~ L 



(2) 



with 7 = 0.26 ± 0.02. This result is consistent with the 
expected DP value 7 = j3/v± (the best known estimate 
2/1 = 1.096854 ± 0.000004 £TJ gives a k 0.252). We 
also estimated r for relatively short times in very large 
substrates (L = 65536), and obtained 



r c (L = oo,t) ~t- 7? , 



(3) 



with r\ = 0.160 ± 0.005. This estimate also supports the 
DP equivalence, which gives 77 = /3/vu, where u\\ is the 
parallel correlation length exponent (best known esti- 
mate u\\ = 1.733847 ±0.000006 JTJ]). 
The interface width, defined as 



W(L,t) 



1/2 



(4) 



obeys dynamic scaling involving exponents of DP and 
KPZ theory (overbars and angular brackets in Eq. 3 de- 
note spatial and configurational averages, respectively). 
In order to understand its behavior below the critical 
point, we first show the results at p c and very large sub- 
strates in Fig. 3. The interface width W increases as 



W ~t,p = p c 



(5) 



as a consequence of the finite fraction of growing columns 
in isolated branches and the increasing fraction of blocked 
columns, which give rise to increasingly large heights' 
differences. 

The evolution of the interface width for p < p c is pre- 
sented in Fig. 4, where we plotted In W as a function of 
the scaling variable x = fe"H , with uu = 1.733847 p9| , in 
substrates with L = 4096. There is a transient region for 
t < t cros ~ e^n , in which W shows the rapid increase 
typical of the critical point (Eq. 5). Notice that t cros is 
the characteristic time of correlations in the DP process. 
At t ~ t cros , W crosses over to a KPZ scaling 



(6) 



with (3k = 1/3 in d = 1. Finite-size effects are responsi- 
ble for the reduced declivities in Fig. 4 when compared to 
the asymptotic forms of Eqs. (5) and (6) (strong finite- 
size effects are typical of ballistic deposition models |Fj] ) . 

For long times, finite-size effects lead to the saturation 
of the interface width. The extrapolation of data for 
several p and L, also considering finite-size effects [ 51] , 
leads to 



W s , 



(7) 



with the KPZ exponent olk = 1/2 and the DP expo- 
nent. e - ' 3 is the typical lateral distance between active 
columns, but appears in Eq. (4) as a vertical scaling 
length, accounting for lateral correlations in the rough- 
ness saturation regime. The divergence of W sa t at p c 
indicates the failure of KPZ scaling at criticality. 

For p > p c , the growth process stops when the whole 
surface is covered with B, for any length L. The heights 
of the blocked deposits attain limitting or saturation val- 
ues with average H s , and the interface widths attain sat- 
uration values W s (W s should not be confused with W sat 
for p < pc, since the former is a property of infinitely large 
static deposits and the latter is related to finite-size ef- 
fects in growing deposits). The time for surface blocking 
is the characteristic time of survival of particles in the 
corresponding CP, consequently H s and W s should be- 
have as the parallel correlation length in the absorbing 
phase: 



(8) 



Eq. (8) is confirmed in Fig. 5, where we show linear fits 
of h\H s and lnW s versus log(— e), with p c — 0.20715 
(the same estimate of the growing phase). From fits with 



2 



different values of p c we obtain v» — 1.75 ± 0.05, which 



is also consistent with DP within error bars 19 1. 

Analogous results were obtained in two-dimensional 
substrates. In Fig. 6a we show logr versus loge, with 
p c = 0.4902, which gives = 0.573 ± 0.020 (Eq. 3). 
In Fig. 6b we show laW versus lnx, x = te u ^, for sev- 
eral values of p, considering v\\ — 1.295 Again it 
shows exponents consistent with DP and the crossover 
from DP to KPZ scaling. At p = p c we obtained Eq. 
(5) with a ~ 0.46, to be compared with the DP value 
a « 0.451 p2] |. The results in d = 2 are less accurate 
due to the limitations in lattice lengths (L < 256), but 
are essential to justify any comparison of our theory with 
experiments. 

The applicability of our model to real growth processes 
is limitted due to the ballistic aggregation conditions, the 
absence of diffusion mechanisms etc. However, if poison- 
ing effects lead to a transition to a blocked phase and if it 
can be interpreted as a transition to an absorbing phase, 
then the robustness of the DP class |20||| suggests this 
type of transition. A possible realization is Si deposi- 
tion in atmospheres with phosphine (PH 3 ), which shows 
a decrease of growth rate with increasing phosphine flux. 
The saturation of phosphorous dangling bonds by hydro- 
gen at the surface was suggested as the main blocking 
mechanism || , but to our knowledge no blocking transi- 
tion was found yet. Another possible application is dia- 
mond CVD in atmospheres with boron, in which the for- 
mation of an amorphous BCN phase blocks the growth 
of the diamond phase for boron to carbon ratios above 
BJC — 0.1 [|J. It seems to be an absorbing transition 
similar to our model and, consequently, is a candidate to 
the DP class. 

Finally, it is important to recall the differences be- 
tween the transition found in our model and the pinning 
transitions by directed percolation of growing interfaces 
in disordered media |j 10 2^Q. In that case the inter- 
face is blocked if the impurity concentration exceeds the 
DP threshold, then infinite surface growth is found in 
the absorbing phase of the impurities system. Conse- 
quently, the critical behavior of geometric quantities such 
as growth rate and interface width are completely differ- 
ent; for instance, Eq. (3) is obeyed with j3 = v\\ |p3| . 
A very different correspondence to DP is also found in 
models with competition between aggregation and des- 
orption that show roughening transitions p"l| ], in which 
the film growth regime parallels the absorbing DP phase. 
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FIG. 1. (a) Examples of deposition attempts in d = 1, 
in which only the configurations of the incident column and 
of neighboring columns are shown. Open squares represent 
particles A, filled squares represent particles B and crossed 
squares represent incident particles (A or B). In processes 
(1), (2) and (3), aggregation occurs at the positions marked 
with a filled circle. In processes (4) and (5) the aggregation 
attempt is rejected. Notice that, in processes (3) and (4), 
lateral aggregation to the right is not possible because the 
neighboring A is not at the top of the column, (b) The equiv- 
alent one-dimensional contact process, in which a top A corre- 
sponds to a particle (empty circles) and a top B corresponds 
to a hole (underlined empty site). The initial configuration 
and the possible final configurations (for the cases of incident 
A or incident B) are shown. 



agencies) . 



FIG. 2. (a) Deposition rate r 
cidence of particles B, in d — 
p c = 0.20715. 



versus probability p of in- 
1. (b) Scaling of r near 
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FIG. 3. Time evolution of the interface width W at the 
critical point p c — 0.20715 in d = 1, for a very large substrate 
(L = 65536). 

FIG. 4. In (W) versus In (a;), with the scaling variable 
x = te"l\ . From below to above, p — 0.15, p = 0.17 and 
p = 0.18 (L = 4096). The regions of critical DP and KPZ 
behaviors are indicated. 



FIG. 5. Saturation height ln(H s ) (squares) and satura- 
tion width ln(W s ) (crosses) versus In— e in d — 1, with 
p c — 0.20715. The solid line is a least squares fit of H s data, 
giving a declivity v\\ m 1.75. 



FIG. 6. (a) Scaling of the deposition rate r near 
p c = 0.4902 in d = 2, giving f3 = 0.573. (b) In (W) versus 
ln(x), with x = te" n , for p = 0.44, p = 0.46 and p = 0.47 
from below to above (L = 256) . 
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